Dynamics in the Metabasin Space of a Lennard-Jones Glass Former: Connectivity and 

Transition Rates 
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Using simulations, we construct the effective dynamics in metabasin space for a Lennard-Jones 
glass-former. Metabasins are identified via a scheme that measures transition rates between inherent 
structures, and generates clusters of inherent structures by drawing in branches that have the largest 
transition rates. The effective dynamics is shown to be Markovian but differs significantly from the 
simplest trap models. We specifically show that retaining information about the connectivity in 
metabasin space is crucial for reproducing the slow dynamics observed in this system. 
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I. INTRODUCTION 

The origin of the dramatic slowing down of dynamics 
in supercooled liquids upon approaching the glass tran- 
sition temperature, has been of great research interest. 
The potential energy landscape (PEL) has proved to be 
an important conceptual tool for analyzing the dynam- 
ics of supercooled liquid [E [3, S II S H ■ In PEL pic- 
ture, the configuration hyper space can be seperated into 
a collection of potential energy valleys, each indcntificd 
with a local minimum or "inherent structure" . There 
have been many recent studies demonstrating that the 
low temperatures dynamics of model supercooled liquids 
is dominated by activated dynamics between "traps" rep- 
resented by metabasins which are clusters of inherent 
structures^, H, Q. In this activated regime, a natural 
mechanism for glassy dynamics is provided by the trap 
model[10]. Extensive numerical simulations of model 
glass formers have shown that the dynamics in metabasin 
space can be mapped on to the original[ll| or an extended 
version of the trap model [12]. The trap model framework 
neglects the connectivity in the metabasin space. An in- 
teresting question to ask is whether a non-trivial connec- 
tivity in the metabasin space, such as the one found in 
the multiple funnel landscape [ll, [13], leads to important 
differences from the trap model predictions regarding dif- 
fusion constants and relaxation timescales. 

An important ingredient in any framework for map- 
ping the dynamics of a supercooled liquid to metabasin 
space is the method for constructing the metabasin. In 
previous work, metabasins were extracted from molecular 
dynamics trajectories projected on to the inherent struc- 
ture spacepjj, where it was observed that there were 
long stretches where the system transitioned back and 
forth between a finite number of inherent structures, and 
these were grouped into metabasins. The usefulness of 
the metabasin concept lies in the enormous simplifica- 
tion in the dynamics that results from projecting the 
supercooled-liquid dynamics on to the metabasin space: 
this projected dynamics is [l2], and therefore, under- 
standing glassy dynamics reduces to the problem of un- 
derstanding the properties of this random walk. Adopt- 



ing the perspective that the aim of the metabasin con- 
struction is to define a space in which the complex dy- 
namics of the supercooled liquid becomes Markovian, 
opens up the possibility of constructing metabasins us- 
ing other algorithms. In this paper, we pursue this line 
of reasoning and use a metabasin construction scheme 
based on a knowledge of the transition rates between dif- 
ferent inherent structures. With the metabasins defined, 
we (a) test whether a trap model with these metabasins 
as the traps provides an adequate description of the ob- 
served dynamics of the supercooled liquid, (b) include 
the connectivity in the metabasin space to extend the 
trap model predictions and show that the correlations 
extracted from this random walk, which takes into ac- 
count the non-trivial connectivity, agrees well with the 
ones measured directly in the molecular dynamics simu- 
lations, and finally (c) we analyze the properties of this 
random walk in order to gain some insight into the slow 
dynamics. 



II. MODEL 

The model system we use to study the dynamics of 
supercooled liquid is a Lennard-Jones binary mixture 
(LJBM) consisting of 67 particles, 53 of type A and 14 
of type B. The simulation box is a cube with period- 
ical boundary condition. All particles have the same 
mass m, and they interact via a Lennard-Jones poten- 
tial v(r) = 4e((cr/r) 12 — (cr/r) 6 ), where r is the distance 
between them. The interaction parameters, depending 
on the types of participating particles, are: caa = 1.0, 
oaa = 1-0, cab = 1-5, (Jab = 0.8, e B b = 0.5 and 
a bb = 0.88[15j. In the following, all quantities will be 
expressed in reduced units, with the unit of length as 
o AAi the unit of energy as €aa, and the unit of time as 
(m<j 2 AA / ^aa) 1 / 2 fl6j |. The number density of particles is 
1.2. To accommodate the box with periodic boundary 
conditions, the potential is shifted and truncated with a 
quadratic cutoff [lfl [l?], [H, EH , which ensures continuity 
of the potential and its gradient at the cutoff radius. 
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The modified potential is 



v c (r) =v(r) — v{r c ) 



(r 2 — r 2 ) ( dv(r) 

-Ae{{a/r) 12 - (a/rf 
+ (6(a/r c ) 12 -3(a/r c ) e )(r/r c ) 2 
-7(a/r c ) 12 +4(<j/r c ) 6 ) 



(1) 



for r < r c and v c (r) = for r > r c . Since the cutoff dis- 
tance should be smaller than half of the box size, a cutoff 
distance r c — 1 .9 is used. This cutoff makes the potential 
shallower than the original and, therefore, the ratio 
between A and B particles is adjusted to minimize the 
chances of crystallization. The velocity form of Verlet 
algorithm is used to solve the Newtonian equation, with 
St = 0.001. The temperature range considered is from 
T = 0.7 to T = 0.48. In this regime the supercooled liq- 
uid slows down significantly yet its equilibrium properties 
can be studied in simulations. Temperature is fixed by 
resampling the velocities from a Boltzmann distribution 
after every 10240 simulation steps. Initial configurations 
are generated by first equilibrating the system at T = 5, 
then at T = 2, and then slowly cooling down to the tem- 
peratures of interest with cooling rate < 3.33 x 10~ 6 [2(| 
(Smaller cooling rates are used at lower temperatures). 
At each of these temperatures, the system is equilibrated 
for a time much longer than (at least an order of magni- 
tude) the estimated a-relaxation time before any data is 
collected for measurement. 



III. INHERENT STRUCTURES 

Inherent structures (IS) are obtained using con- 
jugate gradient minimization techniques during the 
simulation 0, Q. Although, theoretically, inherent struc- 
tures can be labeled by their extact potential energy, 
there is a danger of labeling two different inherent struc- 
tures with the same potential energy since any numerical 
procedures measure the potential energy with hnitc preci- 
sion. To alleviate this problem, we use a pair of energies 
to label each inherent structure, {V,Vbb}i where V is 
the potential energy of local minima, and Vb b is the po- 
tential energy between type B particles. If two different 
inherent structures have the same V, their Vbb are likely 
to be different since they have different arrangement of 
particles. Implementing this procedure is most important 
for the metabasin cosntruction since many thousands of 
IS are generated in the process. 

Given an inherent structure, an ensemble of configu- 
rations in the same valley can be constructed by using 
a restricted Boltzmann sampling that is discussed in the 
next section. With each of these configurations, and ini- 
tial velocities sampled from the Boltzmann distribution, 
the waiting time out of this IS valley can be measured 
using interval-bisection method 0] . For the temperature 
range considered in this paper, it turns out that for many 




FIG. 1: Illustration of the construction of metabasins. Small 
circles denote inherent structures. Solid arrows denote the 
strongest transition branches from corresponding metabasins. 
Each big circle encloses a cluster of inherent structures, i.e., 
a metabasin. Dash lines denote the weaker inter-metabasin 
transitions 



inherent structures, the distribution of waiting time de- 
viates from exponential significantly, thus the transition 
between inherent structures is history-dependent. The 
history dependence indicates that there is no clear sepa- 
ration of scales between the time taken to equilibrate in 
an IS valley and the hopping between different IS. This 
is reasonable given the fact that there can be arbitrar- 
ily small dynamical barriers between the IS, which are 
defined solely based on their property of being a poten- 
tial energy minimum. It is expected that grouping IS 
into metabasins such that IS with frequent transitions 
between them are in the same metabasin would lead to 
larger barriers between metabasins, and, therefore, a sep- 
aration of time scales, and Markovian dynamics. 



IV. CONSTRUCTION OF METABASINS 

Starting from one inherent structure in a hypothetical 
metabasin, the system will more likely go into another 
inherent structure in the same metabasin. To construct 
a metabasin from a randomly choosen inherent structure 
A, we, therefor, need to start the simulation from A for 
many times, and count how many times the system goes 
into each neighboring inherent structure. There must be 
a most frequently visited neighbor, B. Or in another 
word, A has stronger connection to B, than to any of its 
other neighbors. As shown in Figure [TJ A and B get con- 
nected by an arrow. The procedure is repeated starting 
with B and the process of building this cluster contin- 
ues until we find a connection which links two IS that 
are already in the cluster. All the IS in the cluster are 
now assigned to one metabasin. We then check the other 
neighboring IS found in the process of constructing the 
cluster and identify their most frequently visited neigh- 
bor. If these belong to the existing cluster, then these IS 
also get assigned to the metabasin containg the starting 
configuration, A. During the construction of metabasin, 
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tens of thousands of inherent structures are encountered, 
which crowd into a small range of potential energy, and 
as mentioned earlier, it becomes essential to label them 
with both {V,V B b}- 

Our construction is similiar to the hierarchical mas- 
ter equation approach in pll ] . Comparing with MB con- 
struction methods based on trajectories, our approach is 
objective except for the fact that the computational force 
is limited to explore all the details of a MB in the phase 
space. 

The construction of the metabasin, described above, 
relies on the transitions between IS being history inde- 
pendent. It is, therefore, important to obtain an initial 
ensemble of states which are in equilibrium in the valley 
A and use these as starting configurations to measure 
the number of transitions to other IS. All of these start- 
ing configurations should have the same set of values for 
{V, Vbb} but a different set of velocities. We implement 
this restricted Boltzmann sampling by starting from an 
initial configuration Aq, with initial velocities sampled 
from the Boltzmann distribution, and run the simula- 
tion for n MD steps. If the final state Ai still belongs 
to the same inherent structure, then we accept it as a 
member of the ensemble. If not, we count Aq again as 
a member of the ensemble. This process is repeated un- 
til we have a large enough ensemble of initial states to 
measure the transitions from. In practice, n is chosen 
to be 1/4 of the estimated waiting time of the given IS, 
and in order to reduce correlation, the above process is 
repeated 8 times before the resulting state is accepted as 
initial state. The number of initial states is at least 40, 
and more often chosen to be greater for better estimation 
of branching ratios. 

The temperature chosen to construct metabasin is 
T = 0.52, a temperature where IS arc well defined and 
the dynamics is not prohibitively slow for metabasin con- 
struction. The number of initial inherent structures, ran- 
domly chosen from MD trajectories is 140 and each of 
these constitute the starting point for metabasin con- 
struction according to the algorithm described above s. 

The potential energy, E, of a metabasin is defined to 
be the energy of the most probable inherent structure in 
it. 



V. MARKOVIAN PROPERTY 

A hallmark of Markovian processes is that the distribu- 
tion of waiting times is exponential [22| . Since the motiva- 
tion behind the construction of the metabasins is obtain- 
ing a Markovian model of the dynamics of supercooled 
liquids, the first task is to check whether our metabasin 
construction yields a space in which the dynamics is 
Markovian. To measure the waiting times, an ensemble 
of initial states is constructed through a restricted Boltz- 
mann sampling in each metabasin and at the temper- 
ature of interest. A constant temperature trajectory is 
then started from each initial state, and interval-bisection 
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FIG. 2: Effective activation energy E act of metabasins, plot- 
ted against their energy per particle, E/M, where M is the 
number of particles. The activation energies are calculated 
according to Eq. 



method [9( is used to determine the inherent structure se- 
quence of the trajecotry. If the inherent structures of two 
succesive quenchs are not in the original metbasin, the 
system is considered to have made a transition to a new 
metabasin, and the waiting time is recorded. Since the 
conjugate gradient method used to minimize the poten- 
tial energy occasionally leads to a wrong inherent struc- 
ture, two succesive quenchs are used to signal a transi- 
tion. For all the metabasins constructed according to our 
algorithm, the waiting time distributions are founded to 
be exponential and, therefore, the dynamics in this space 
is indeed Markovian. The mean waiting time Tj, for a 
metabasin i with energy Ei is measured in the temper- 
ature range 0.44 < T < 0.70. This mean waiting time 
is found to be above the ballistic time region, indicating 
the separation of time scales that we would like to see for 
metabasins. The mean waiting time is characteristic of 
an activated process with: 



(2) 



The effective activation energy E act is not strongly cor- 
related with the energy E although E act tends to be 
higher for lower lying metabasins, as shown in Figure [5] 
The time scale, t(T) demonstrates a stronger correlation 
with a marked increase for deeper metabasins, as shown 
in Figure [3] Since the dynamics in metabasin space is 
Markovian, the evolution of the probability distribution 
of metabasins, P,(t), is described by: 



1 



Ti{T) 



(3) 



where 7$ is the mean waiting time of metabasin i at tem- 
perature T, and Wij is the transition rate from metabasin 
j to metabasin i. Envisioning the metabasins as traps of 
depth Ef ct , the simplest trap model is one that ignores 
the connectivity and characterizes the trap space by n 
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and a density of states [l0|: 
±P(t) = - — 



10" 



Pi(t) + u(t)p(E?*) 



(4) 



where w(t) is a normalization factor and p(E act ) is the 
density of traps with a depth E act . In most mappings of 
metabasin dynamics, in Lennard- Jones systems, to trap 
models the energy of the metabasin, E has been consid- 
ered as the depth of the trap. As discussed earlier, we 
find that the activation energy and the nominal energy 
of a metabasin are not strongly correlated. We, there- 
fore, choose to construct an effective trap model that 
describes dynamics in an energy landscape, by averaging 
Ti over metabasins with energies in a small range around 
a value E: 



dt 



P(E,t) 



t(E,T) 



p(E,T)u(t) 



(5) 



where 



P(E,t) = '£5(E i - E)Pi(t) 

1 J2 i 6(E i -E)P i<eq (T)/r i (T) 



r(E,T) 
P(E,T) 



Y, i S(E i -E)P i , eq (T) 
P eq (E,T)/r(E,T) 



(6) 



with Pi. e q(T) oc exp(— (3Ei), the equilibrium probability 
of finding metabasin i, and P eq (E, T) the equilibrium dis- 
tribution of metabasin energy, which can be measured in 
the simulation. As shown below, this effective trap model 
differs from the simple trap model 10] in that p{E, T) de- 
pends on temperature; reflecting the fact that Ei is not 
the trap depth of metabasin i. 

The metabasin energy distribution P eq (1?, T) is ob- 
tained from simulations run at T = 0.6, 0.52 and 0.48, 
and is found to be well described by a Gaussian over 
this range 0, H, [ll], [23|. We have extracted the density 
of states from these distribution and find that fl(E) oc 
P eq (E,T) exp(i? /T) is, indeed, independent of temper- 
ature as show in Figure We use this observation to 
construct P eq (E, T) at temperatures other than the ones 
where it was measured explicitly, and combine this infor- 
mation with t(E,T) to obtain p(E,T). This function, 
shown in Figure 21 depends on temperature and is dif- 
ferent from fl(E). 

The effective trap model based on the deduced forms 
of p(E,T), t(E,T) and u>(t)), provides a description of 
the activated dynamics near the glass transition. It is, 
however, a meanfield model that ignores the connectivity 
in the metabasin space. In order to evaluate the effects 
of the connectivity on the dynamics, we have calculated 
various correlation functions using this trap model and 
compared it to the results of the direct, molecular dy- 
namics simulations. 
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FIG. 3: Waiting time of metabasins measured at T = 0.6 and 
0.48, plotted against their metabasin energy per particle. The 
solid line is t(E,T) calculated according to Eq. © 
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FIG. 4: p(E, T) at T = 0.6 and 0.48, calculated according to 
Eq. © 



VI. TESTING THE TRAP MODEL 

In this section, we compare the predictions of the effec- 
tive trap model to actual simulation results. In the trap 
model, correlation functions are calculated by averaging 
over the sampled metabasins according to: 

{A) eq = Y d Pi,e q {T)A i 

J J2 t ^e q (T)S(E-Ei) 
^E, T) \%^ (7) 

Here A is the physical observable of interest. The first 
summation, J^. is over all metabasins, while ^ Q is over 
the metabasins, sampled according to the Boltzman dis- 
tribution in the simulations. As mentioned earlier, a total 
of 140 metabasins were sampled. In numerical calcula- 
tion, the 8 function in ^ , is replaced by a Gaussian with 
a narrow width, 2a 2 — 0.01 2 M 2 , with M = 67 being the 
number of particles. 
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FIG. 5: Plot of the density of states of metabasins, Q(E), 
shifted vertically to make lines overlap. 

A. Estimate of the alpha relaxation time 

The a relaxation time r Q was measured in the sim- 
ulation through its usual definition: F s (qo,T a ) = 1/e , 
where F s (q,t) is the incoherent scattering function [24|. 
In this paper, F s (q,t) is measured for type A parti- 
cles and qo = 7.251 [151 ]. The a relaxation time can 
also be estimated from the distribution of single particle 
displacements [25j, which shows marked non-Gaussianity 
at timescales of the order of r Q . The non-Gaussianity 
is related to caging and a caging time scale, r s can be 
extracted by measuring the time at which the probabil- 
ity of a type A particle having a translation less than 
l/<?o diminishes to 1/e. As shown in Table[H in the tem- 
perature range examined, t s and r Q are proportional to 
each other. The caged motion of particles has also been 
related to hopping between metabasins Q. 

In the trap model framework, the a relaxation process 
corresponds to hopping between metabasins 0, lll(, and 
in this regime, the correlation function C(t) can be con- 
structed by assuming that the correlation is unity when 
the system is in the same metbasin and drops to zero as 
it leaves the metabasin[lj|- For our effective trap model, 
this approximation leads to: 

C(t) = (exp(-t/ri))e« (8) 

The expectation value is calculated according to Eq. |(7J). 
In Figure El C(t) is plotted along with F s (q,t). As ex- 
pected the shapes of the two functions are different since 
C(t) has no information about the short time dynamics. 
For C(t) to be a useful tool for understanding the glass 
transition, however, the change in time scale of the long- 
time dynamics should closely resemble that of r a . Within 
the trap model, a measure of the a relaxation time is pro- 
vided by the relaxation time r c defined by C(r c ) = 1/e, 
and measured using Eq. ([5]) . The values of r c along with 
the values of r a are listed in Tablc|T]for T = 0.6, 0.52, and 
0.48 and plotted in Figure [7] It is clear that r c is much 
smaller than r a for all temperatures listed, and more im- 
portantly, r a increases significantly faster than r c . These 



FIG. 6: Comparison of C(t) obtained from the effective trap 
model at T = 0.6 and 0.48 and the incoherent scattering func- 
tion, F 3 (q,t) measured in simulations. Also shown is F' s {q,t), 
the incoherent scattering function measured from the inher- 
ent structure coordinates. Time is measured in scaled units 
(c/ text). 



TABLE I: Comparison of estimates of the a relaxation time, 
r a , t s , and T c . 



T 


Ta(T) 


Ts{T) 


Te(T) 


Tc{T) 
To (0.6) 




rc(T) 
t c (0.6) 


0.9 


2.2 


6.7 




0.11 


0.11 




0.6 


20 


62 


2.9 


1 


1 


1 


0.52 


91 


2.7 x 10 2 


1.1 


4.5 


4.4 


3.8 


0.48 


3.7 x 10 2 


1.1 x 10 3 


32 


19 


18 


11 


0.46 


8.4 x 10 2 


2.5 x 10 3 




42 


40 





results indicate that the trap model, as defined up to now, 
does not capture all of the physical processes leading to 
the slow dynamics. 
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FIG. 7: Temperature dependence of T a and r c , measured in 
scaled units. The solid line is the fit of r a to the Vogel- 
Tamman-Fulcher form. 
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FIG. 8: MSD of type A particles plotted as a function of time 
t measured in scaled units. The lines correspond to T = 1.5, 
0.9, 0.7, 0.6, 0.56, 0.52, 0.48, and 0.46, top to bottom at t > 1. 



B. Trap model prediction of the mean squared 
displacement (MSD) 

The MSD (5r 2 (t)) of type A particles is measured di- 
rectly from the simulations, and as shown in Figure [5J 
there is caging. The diffusion constant can be measure 
by looking at the long-time behavior: D = (Sr 2 (t))/6t 
at large (5r 2 (t)}. The temperature dependence of D is 
shown in Figure [TOl 

The MSD of particles in space can be mapped on to 
a random walk between metabasins 0, [2(| . Consider a 
random walk going through a sequence of N metabasins, 
ot\ — > CK2 — > ••■ — > &N+1, where a n is the label of 
metabasin. The metabasin coordinate sequence is labeled 
as £i — > £2 — > ••• — > £zv+i- where £ n is a 3M — component 
vector representing the coordinates of M particles in 
metabasin n. The diffusion constant, Dmb, of this ran- 
dom walk can be measured from the large N limit: 



6MD MB 



JV+l 



£i) 2 } ? 



N(r n } 



(9) 



where (r n ) ru) is the average waiting time obtained 
by averaging over all random walk trajectories in the 
metabasin space. For each metabasin a n , the probabil- 
ity of being at metabasin i is proportional to Pi^ql^i- 
Thus (T n )rw can be calculated from (l/ri) e g, which is an 
average over all metabasins weighted with equilibrium 
probabilities, 



Pi,e 



1 



(10) 



<■</ 



Assuming that at each metabasin transition the displace- 
ment in £-space is vA, the diffusion constant Dmb is, 
therefore, Dmb — (A/6M)(1/Vi) e g. As shown in Fig. 
flUl Dmb rides above the D measured from the MSD 
of particles, again, demonstrating that the trap model 
picture, which ignores all connectivity in the metabasin 
space, underestimates the slowing down of the dynamics. 



VII. EFFECTS OF CONNECTIVITY IN 
METABASIN SPACE 



One feature of the metabasin space that has been ig- 
nored in the trap model is that the Markovian dynam- 
ics describes a random walk on a graph with non-trivial 
connectivity. The calculation of the diffusion constant, 
Dmb can be improved upon by incorporating aspects 
of the connectivity that are captured by the transition 
rates Wij. Defining, 8£ n — Cn+i — £n, the mean squared 
displacement in £-space can be written as: 

N 



((6v+l - £l) )rw = 5£ n ) ) rw 

n=l 

N N 

n— 1 n—2 
N-l N 



rn — 2 n— TTi+1 

N(A + 2(5£ n -i ■ O™) 



(11) 



Here we ignored (<5£„_ m ■ <5£„ 
expression of Dmb becomes 



for m > 1. Then the 



6MD MB « (l/r<)e,(A + 2(^„_ 1 • 8& 



(12) 



To calculate (<5£„_i • 8^ n ) rw , one first picks the 
metabasin a n at time step n to be a specific one, say 
i, and averages over all possible metabasins j at step 
n — 1, and over all possible metabasins j' at step n + 1, 
followed by averaging over all possible metabasins i at 
step n. For a particular a n = i, the probability of having 
a metabasin j' at 71 + 1 is Wj'j/Efc w ki — w j'i T i- Simi- 
larly, the probability of having a metabasin j at step n — 1 
must be 10,-jTj. Thus 



(<&-i>< ■ (<C 



j'i) 



(13) 



with 5£ij = £j — Two facts should be noticed, (1) 
the dimension of £, i.e., the number of independent co- 
ordinates, d = 3M — 3 = 198, is very high, and (2) For 
many metabasins the connectivity is sparse, with strong 
connections only to one other metabasin. An example is 
provided by the metabasins A and B in Table HT1 Because 
of (1), any pair S^ij and 5£j*i with j' ^ j can be consid- 
ered as two random vectors in a high dimensional space, 
and hence, S^ij -S^jn f=a 0. The exception is of course when 
j = j' , i.e., when the system hops back and forth between 
two metabasins i and j, and <5£ij i = — (S^ji) 2 = — A. 
So we have 



(££„-i • <S£ n )i ~ - y^(wj»Tj) 2 A 
3 

= -P 2 A 



(14) 
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TABLE II: The branching ratio, 
motabasins at T = 0.52 



in of four different 



i 


A 


B 


C D 


UJj 1 1 % 


0.87 


0.6f 


O.ff 0.f3 




0.05 


0.24 


0.06 0.07 




0.02 


0.06 


0.05 0.05 




0.007 


O.Of 


0.04 0.04 




0.006 


O.Of 


0.04 0.03 




0.003 


O.Of 


0.04 0.03 




0.76 


0.43 


0.03 0.04 



where pf = J2j( w ji T i) 2 ^ the probability of return to 
the previous metabasin. For those metabasins that have 
many neighbors, each with similiar (and low) branching 
ratio Wj'iTi, such as metabasin C and D in Table HT1 the 
return probability is small. 

Averaging Eq. (fl~4"]) over all possible metabasins i, one 
obtains, 



{ l / T i)eq 

and the expression for the diffusion constant becomes 



(15) 



D 



MB 



A 

6M 



((1-2^)M) 



eq 



(16) 



The branching ratios WjiTi are measured for all the 
sampled metabasins, and pfs are calculated. As shown by 
the examples in Figure^ the return probability pf is rela- 
tively low at high temperatures. But, as the temperature 
decreases, pf increases significantly for some metabasins. 
This is because some of the transition branches are sup- 
pressed as temperature decreases. According to Eq. ([IB"]) , 
high pf makes the diffusion constant smaller. Increasing 
values of pf is , therefore, another source of the slow dy- 
namics, and one that is not captured in the trap model. 

In Figure \TU[ the prediction of Eq. (fTTj)) is compared 
to measured D(T). It shows that incorporating the re- 
turn probability pf has a significant effect on Dmb and 
improves the agreement with the measured D. In Q, 
Doliwa and Heuer were able to predict the diffusion con- 
stant without considering the return probability, due 
to the fact that the metabasins were constructed from 
trajectories, and the inherent structures that were con- 
nected by quick hops were grouped together into the same 
metabasin. The dynamics in the space of MB defined in 
this manner satisfies more stringent requirements [27J, in- 
ducing pf = 0. In this work, metabasins were instead 
constructed from the transition branch ratio of inher- 
ent structures, without any consideration of the absolute 
values of the barrier heights between metabasins. Our 
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FIG. 9: Comparison of pf(T), the return probabilities of the 
four metabasins listed in Table ITTl 
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FIG. 10: Comparison of the temperature dependence of the 
measured diffusion constant, D, of type A particles, with the 
predictions of the trap model, and the trap model modified to 
include a finite return probability. The relationship between 
the diffusion constant and n in Eq. [16] has an undetermined 
constant, which was chosen to be 1/340 in this figure. 



construction gives rise to a graph which is much more in- 
homogeneous and the mean-field assumption of the trap 
model is not a good approximation. The advantage of 
this method of construction is that it gives a better idea 
of the connectivity. To emphasize, the effective dynam- 
ics in the metabasin space is certainly Markovian but the 
random walk lives on a graph with non-trivial connectiv- 
ity. 



VIII. SUMMARY 

In this work, we have presented a technique for con- 
structing metabasins using branching ratios of transi- 
tions between inherent structures. The dynamics in the 
space of these metabasins is Markovian, and can be de- 
scribed as an activated process in an energy landscape. 
A mean field trap model that ignores the connectivity 
of the metabasins, however, underestimates the slowing 
down of the dynamics with decreasing temperature. We 
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show that incorporating the connectivity through a set 
of return probabilities leads to quantitative agreement 
between the predicted and measured temperature depen- 
dence of the diffusion constant. This result demonstrates 
that connectivity is an essential element of the effec- 
tive dynamics in the metabasin space. At temperatures 
higher than ~ 0.7, the mean waiting time of metabasins 
approachs the ballistic time scale, and mapping the dy- 
namics to metabasin space loses its utility. At tempera- 
tures lower than ~ 0.45, the return probabilities become 
so large that higher order terms in Eci. tTTTj) need to be 
included and the effective dynamics loses its simplicity. 
The mapping to metabasin space is, therefore, most use- 
ful in the range 0.45 < T < 0.7 At the high end of this 
temperature range, the connectivity is reasonably homo- 
geneous, and the return probabilities are low. In this 
regime, effective trap models capture the essential fea- 
tures of the slow dynamics. As the system is cooled, some 
of the connections are suppressed more strongly then the 
others. The connectivity becomes sparse, resulting in 
higher return probabilities, and break down of the as- 
sumption, underlying the trap model, that once the sys- 
tem escapes fro a trap it chooses one from a fixed distri- 
bution of trap depths that does not have any connectivity 
information. It is clear that the difference between con- 



nections, i.e., the inhomogeniety of connectivity causes 
the connectivity to be sparse at low temperatures. At 
the same time, we find no obvious correlation between 
return probability and metabasin energy, and therefore, 
metabasins of different energies seems to have similiar 
connectivity environment. 

Instead of metabasins, J. Kim and T. Keyes [28J have 
considered the strong return dynamics between inherent 
structures of a supercooled CS2 system, and postulated 
that the return dynamics could be coarse-grained to the 
motion within metabasins. In this paper, however, we 
showed that even after coarse-graining the strongly con- 
nected inherent structures into metabasins at a temper- 
ature in the middle of the range of interest to the slow 
dynamics, the connectivity between metabasins plays a 
significant rule, especially at low temperatures. Return 
dynamics appears to be less important for some small 
systems [29f. In general, however, if the multiple funnel 
landscape[l3l [l4| is an adequate representation of the en- 
ergy minima, then any definition of metabasins based on 
a single temperature will have to include the influence of 
return dynamics at low enough temperatures. 

We acknowledge helpful discussions with David Reich- 
man, especially during the early stages of this work. This 
work was supported by NSF DMR-0549762. 
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